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This  report  presents  an  analytical  method  for  calculating 

residual  stresses  in  joints  riveted  dynamically,  results  of  the 
! 

Moire  experimental  method  for  measuring  the  residual  strain 
distribution  around  a rivet  hole,  and  describes  an  experimental 
method  for  measuring  the  radial  velocity  of  an  expanding  rivet. 

The  advantages  of  compressive  residual  stresses  in  dynamically 
riveted  joints  are  demonstrated  through  constant  amplitude  fatigue 
testing. 

N. 

The! paXticular  cases  in  which  residual  stresses  were  inves- 
tigated include  aluminum  7075"T6  and  2024-T3  sheet  with  in. 
diameter  holes*,  riveted  with  A"286  rivets  installed  with  the 
Stress  Wave  Riveter.  Stress  wave  riveting  expands  the  rivet 
radially  under  high  acceleration  to  achieve  favorable  stress 
distributions.  Analytical  predictions  of  these  stress  dis- 
tributions are  compared  with  experimental  results  obtained  using 

I 

the  Moire  photoelastic  method.  To  demonstrate  the  benefits  of 
compressive  residual  stresses  in  riveted  joints,  fatigue  specimens 
made  of  2024-T8i  aluminum  were  used.  The  specimens  were  tested 
at  constant  amplitude  load  cycling.  In  the  first  series  of 
tests,  the  fatigue  life  of  specimens  with  residual  stresses 
was  approximately  ten  times  higher  than  the  fatigue  life  of 
specimens  with  open  holes.  In  the  second  series,  precracked 
specimens  with  a 0.05  in.  crack  in  each  hole  were  riveted  with 
the  Stress  Wave  Riveter.  The  cracks  were  oriented  so  that  the 
fatigue  loading  was  normal  to  the  crack  path.  The  fatigue  life 
of  these  specimens  was  also  approximately  ten  times  higher  than 
the  fatigue  life  of  those  specimens  with  precracked  holes  without 
rivets.  igy 
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1.  INTRODUCTION 

Fatigue  life  of  structures  Is  usually  determined  by  the 
growth  of  cracks  originating  at  points  of  local  stress  concen- 
tration. The  holes  provided  for  fasteners  in  riveted  joints  are 
known  to  be  the  most  probable  locations  for  crack  initiation 
due  to  a)  high  frequency  of  defects  in  the  hole  machined  surface 
and  b)  high  concentration  of  in-service  tensile  stresses  near  the 
hole  svirfaoe.  Substantial  and  significant  benefits  in  design 
life  and  structural  weight  can  be  achieved  by  the  use  of  special 
techniques  to  reduce  the  levels  of  maximum  tensile  str€iss  near  a 
hole  by  introducing  residual  compressive  stresses  that  prevent 
or  inhibit  the  growth  of  cracks  under  cyclic  loading. 

It  is  known  that  plastically  deforming  (i.e.  cold-working) 
the  hole  material  Induces  compressive  residual  stresses.  Dietrich 
and  Potter  (Ref.  1)  used  x-rays  to  show  the  existence  of  residual 
compressive  stresses  that  resulted  from  the  use  of  a mandrel. 

Chang  (Ref.  2)  showed  that  such  residual  stresses  do  inhibit 
crack  growth  and  improve  fatigue  life  in  riveted  joints.  This 
report  discusses  the  prediction  and  experimental  verification 
of  the  residual  stress  distributions  that  are  produced  by  dynamic 
installation  of  rivets. 

An  Important  distinction  must  be  made  between  the  following 
methods  of  fatigue-resistant  fastening: 

a)  Precision  fasteners,  or  threaded  pins,  driven  into 
undersized  holes  (i.e..  Taper  Loks,  High  Tigues , etc.) 

b)  Fasteners  that  take  advantage  of  residual  stresses 
induced  in  the  hole  through  burnishing,  peening, 
or  forming  using  a mandrel  or  a hard  steel  ball. 

The  fastener  is  of  secondary  Importance  in  these 
methods. 
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c)  Rivet' type  methods  which  utilize  a rivet  that  is 

usually  subjected  to  higher  loads  than  those 
necessary  to  fill  the  hole,  thereby  producing 
residual  compressive  strei^ses  in  the  surrounding 
metal. 

Precision  fasteners  (a)  require  precision  holes  (typical 
tolerances  of  >-  0.001  in.),  and  after  they  are  installed  there 
is  little  unloading  that  could  allow  the  material  that  expanded 
plastically  to  relax  with  significant  compressive  residual  stresses. 
Fastener,  hole  preparation  and  quality  control  cof.s  are  usually 
quite  high  with  these  methods.  Residual  stresses  are  produced  in 
(b) , where  a section  of  the  hole  is  put  in  a plastic  state  so 
that  the  surrounding  elastic  material  squeezes  that  section  as  it 
does  in  the  "autofrettage"  treatment  of  thick  "wall  cylinders 
(Ref.  3).  This  class  of  processes  is  labor-intensive,  and  has 
not  yet  been  successfully  automated. 

This  report  is  also  concerned  with  the  computation  of 
residual  stresses  and  strains  by  dynamic  methods  of  the  type 
described  in  (c),  which  offers  very  inexpensive  fastening.  Of 
the  types  of  riveting  processes  included  in  (c)  we  shall  consider 
only  dynamic  riveting  where  the  rivet  expands  with  sizeable 
velocity  and  large  radial  accelerations.  One  such  process  is 
implemented  by  the  Stress  Wave  Riveter  (Ref.  4). 

There  are  various  methods  available  for  calculating  residual 
stresses  (Refs.  5,  6).  Most  of  them  are  Implemented  by  computer 
prograras  using  finite  element  schemes.  The  difficulty  with 
these  programs  lies  in  their  inability  to  do  numerical  examples 
with  various  cases  without  extensive  involvement  of  the  user 
toward  understanding  the  program  Itself.  In  this  work  explicit 
equations  are  used  that  provide  a well-defined  series  of  equations 
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that  can  be  solved  on  a programnable  desk  calculator.  The  user 
can  readily  understand  the  cause  and  effect  relationship  in 
every  new  case  without  the  conxnitnaent  inherent  in  large  computer 


2.  DERIVATION  OF  EQUATIONS 


PROBLEM  DESCRIPTION 

In  riveting  with  the  Stress  Wave  tool,  a rivet  is  placed  in 
a hole  with  the  bucking  bar  on  one  side  and  the  Stress  Wave  tool 
on  the  other  (see  Fig.  la).  At  first  the  stress  wave  overcomes 
the  elastic  limit  of  the  rivet  material  and  sets  it  in  plastic 
flow  (Ref.  4).  Subsequently,  as  the  massive  head  of  the  Stress 

Wave  tool  (see  Fig.  lb)  moves  axially  against  the  rivet,  the 

rivet  expands  radially  out  toward  the  surface  of  the  hole.  The 
analysis  in  this  report  is  concerned  with  the  stress  and  strain 
distribution  in  the  material  around  the  hole,  from  the  time 
(t  ■ 0)  when  the  rivet  impacts  the  hole  surface  to  the  time 
(t  ■ t^)  when  the  dynamic  process  of  riveting  terminates.  In 
formulating  the  problem  of  stress  distribution  we  have  hoop 
stress  (o'g),  hoop  strain  (e^),  radial  stress  (a^)  and  radial 
strain  (e^.);  also,  (the  radial  position  prior  to  the 

riveting  disturbance)  and  the  time,  (t).  The  following  two 
boundary  conditions  are  used;  a)  the  velocity  history  at  r^  - a 
(see  Fig.  Ic)  and  b)  the  condition  of  zero  stress  and  strain  at 

some  radius  far  away  from  the  hole  r^  ■ b. 

The  velocity  history  at  r^  « a is  found  experimentally 
(see  Fig.  2).  It  has  three  parts;  a)  linear  loading  from  t « 0 
to  t - t^,  b)  constant  velocity  from  t ■ t^  to  t ■ t2,  and 
c)  unloading  from  t - t2  to  t - t^.  Figure  2 shows  the  entire 
riveting  sequence  by  plotting  strain  vs.  time.  Reference  4 
describes  the  process  in  complete  detail. 
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B)  EQUIVALENT  SCHEMATIC  OF  STRESS  WAVE  RIVETING 


C)  VELOCITY  HISTORY  AT  = a + 


Fig.  1 Corcept  of  Stress  Wave  Riveting 
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Fig.  2 Strain  Gauge  Record  of  1/4  in.  A-286  Rivet 


GENERAL  RELATIONS 


The  relations  used  In  clastoplastlc  problems  of  cylinders 
under  Internal  pressure  are  suitable  for  our  application  with 
rivets.  In  polar  coordinates,  assuming  axlsymmetrlc  conditions, 
these  relations  are  as  follows  (Ref.  7). 


Equlllbrltim  Equation  In  Radial  Direction 


au  , , °r  ~ 

at  ar  r 

o 


o 


0 


(1) 


Compatibility  Equation 


^r 

o 


0 


(2) 


Incompressibility  Equation 


(3) 


Strain  Displacement  Relations  (See  Appendix  A) 


'r  * ^r 


Stress-Strain  Relations 
1 


r E ^ r 


"0  "i  ^^0  ■ ^<^r  + "0 


"z  " i t'^z  ■ ^<^r  + ^0>]  ' <^0  + "r^ 


(4) 

(5) 

(6) 
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For  cases  of  plane  stress  a - 0 and  for  plane  strain  e • 0. 

z z 

In  both  cases  the  shear  stresses  and  strains  are  zero.  The 
plane  strain  condition  will  be  used  In  this  analysis. 

In  our  analysis  we  find  it  convenient  to  use  the  mean 
normal  stress  defined  as  (Ref.  5) 

S - -I  (Or  + ag  + 

where  Oq,  and  are  the  principle  stresses. 

The  principle  components  of  the  stress  deviator  tensor  are 
as  follows 


(7) 


In  the  analysis  we  use  the  deformation  theory  (Ref.  7,  p.  120) 
(in  contrast  to  incremental  or  flow  theories)  in  which  the  plastic 
strains  are  functions  of  the  current  state  of  stress  and  are 
independent  of  the  history  of  loading.  We  are  permitted  to  use 
the  deformation  theory  because  of  the  axial  symmetry  and  radial 
loading  of  the  problems  encountered  in  riveting  and  because  we 
consider  separately  the  loading  and  unloading  history  of  the 
problem. 

The  mean  normal  strain  is  given  by 

« “ 3 ('r  + "e  + <8) 


such  that  Hooke's  law  can  be  expressed  by: 

S * 3k  • e , where  k is  the  bulk  modulus 
The  principal  strain  components  are  given  by 


(see  Appendix  A) 
E 

" 3(l-2^i) 


(9) 
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e and  e 
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Under  the  incompressibility  assumption,  however,  e • 0 
thus 


^r  ' 


€.  and  e 

a Z 


The  following  expressions  of  the  generalized  Hooke's  law 
can  be  used 


9u 


(1+u) 


k‘e  + 


(1-hi)  r 


and 


9u 


(14U) 


k*e  + 


X 


<l+u) 


9u 


k*e  + 


"0  (1-Hi)  (l+u)  "e 

We  observe  that  incompressibility  does  not  permit  us  to  solve  for 

S in  the  equation  3k* e • S,  since,  when  e o,  we  have  k oo. 

From  Hooke's  law,  however,  we  know  that  a + + a • 9k*e 

r “ z * 

which  leads  us  to  conclude  that  + S.  + S - 0.  Furthermore 
from  the  definition  of  S^,  Sg  and  and  Hooke's  law  we  find 

Sr  - e' 


'e 


E e. 


(10) 


and 

where 


E e 


! 

E • (see  Appendix  A). 

Considering  the  case  of  plane  strain  where  e » 

z 

all  times,  we  have 


0 and  S » 0 at 
z 


and 


"r  S 


0 

0 


(10a) 
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CRITERIA  FOR  YIELDING 


In  order  to  extend  our  analysis  to  problems  of  riveting  with 
plastic  deformations  we  must  set  up  criteria  for  yielding.  It 
is  generally  accepted  that  hydrostatic  pressure  does  not  produce 
appreciable  plastic  deformation.  Instead,  it  is  usually  assumed 
that  plastic  deformation  depends  only  on  the  stress  deviation. 

Using  the  well  known  Von  Mises  criterion  (Ref.  7)  for 
yielding,  we  have 


[(s^  - + (Sg  - + (S^  - - K (11) 

where  K is  the  yield  stress  in  pure  shear.  It  can  be  shown 
(Ref.  5),  that  K « where  is  the  yield  stress  in  the 

uniaxial  tension  loading. 

Thus,  with  S^,  • 0 

+s^  Sg  + (12) 

APPLICATIONS  TO  DYNAMIC  RIVETING 

Combining  Eqs.  (1)  and  (3)  we  obtain  a system  of  linear 
first  order  partial  differential  equations  that  describe  the 
initial  (elastic)  motion  of  the  rivet. 


where 
and 

This  equation  is  obtained  by  substituting  the  strain  conditions 

into  the  incompressibility  conditions.  Substituting  for  a and 

r o 

from  Eq,  (7),  and  differentiating  the  second  equation  with  respect 
to  time  we  obtain 


„ ou  , r 

P si 


u 


hr 
" 5t 

^ + -t  . 0 

or  r 
o o 


11 


and 


(Sr  - hi 

r 

o 


0 


(13) 


The  above  two  partial  differential  equations  can  be  reduced 
to  one  partial  differential  equation,  with  one  dependent  variable 
S and  two  Independent  variables  r^  and  t,  by  employing  Eq.  (3) 
to  obtain  and  Eq.  (10)  to  obtain  S^.  and  in  teimns  of  S. 

The  velocity  function  u(t)  at  r^  « a,  and  the  fact  that  S » 0 
at  all  times  at  r^  • b (b  » a) , are  the  two  boundary  conditions 
used  in  the  solution.  Details  are  given  in  Appendix  B. 


The  elastic  solution  for  the  stresses  and  strains  is  as 
follows : 


- pna£n(“) 
o 


E mt  a 


Oq  • pma£n(^)  + 

o o 

mt^a 

^r  " ■ 2r  2 
o 

and 

mt^a 
" 2r  2 


RELATIONS  UNDER  PLASTIC  DEFORMATION 


(14) 


Plastic  deformation  begins  at  r^  - a,  whenever  the  radial 
pressure  Increases  to  the  point  that  S^.  --K  and  Sq  » +K  such  that 
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Von  Mises'  criterion  for  yielding  applies;  that  is,  S S_  S,,  + 
2 2 r r o 

Sg  » K . Substituting  »-K  at  r^  » a in  the  first  Eq.  of  (14) 

we  can  find  the  time  t^  that  yielding  will  commence. 

The  stress  and  strain  relations  in  the  plastic  region  are 
derived  as  in  the  elastic  case.  Since  the  incompressibility 
relation  applies  to  both  elastic  and  plastic  cases,  the  velocity 
relationship  is  the  same,  that  is, 


Considering  materials  that  are  perfectly  plastic,  we  have 


CT^  • S-K,  Oq  • S+K  and  cTj.  " cTq  • 2K 
Substituting  in  Eq.  (1)  we  obtain 

.ana  . JiL  . 2ii  . 0 

r ^ Sr  r • ° 

o o o 

Integrating  we  obtain 


(16) 


S • (2K  - pma)  ^n  r^  + f(t) 

At  the  elastic-plastic  interface  (r^  - R)  both  the  plastic 
and  elastic  solutions  hold. 


Hence , 


- S - K - (2K  - pma)  £n  R + f (t)  - K 


where , R ■ 
interface 


maEt 

2K 


- fma^n  , 

^ 2R^ 

, is  the  radius  of  the  elastic  plastic 


from  which  we  find. 


t 2 

f(t)  - K + pmain  (“)  - ^ - (2K  - pma)  £n  R (17) 

^ 2R^ 

Thus,  the  stresses  in  the  region  a r^  ^ R are  given  by 
aj  - (pma  - 2K)  in 

jg  • (pma  - 2K)  in 

o 

The  stresses  in  the  region  R s:  r^  b are  given  by  Eq.  (14). 

REGION  OF  CONSTANT  VELOCITY 

In  the  two  previous  sub-sections,  the  elastic  and 
plastic  solutions  of  the  stress  and  strain  distributions 
were  found  for  the  loading  phase  where  the  velocity  of  the  hole 
surface  was  a linear  function  of  time  (l.e.,  u » mt).  In  this 
section,  the  solution  is  extended  to  the  region  where  the 
velocity  at  r^  • a remains  constant  (i.e.  , u ■ 

The  solutions  are  derived  in  Appendix  C.  The  elastic 
solutions  (i.e.,  r^  ^ Rj^)  are  given  by 


' 2 

.bv  E rot  a 
+ pnain  (*)  - ^ 


2R‘ 


> (18) 


+ 2K  + proa£n  (J)  - ^ 

^ 2R^  ^ 


14 


and 


mat 


t 

^ i* 


(19) 

(cent) 


The  plastic  solution,  (i.e.  , a ^ ^ Rj^)  are  given  by 


, E mat 

jP  - (xaaZn  (~)  - ^ 

^ ^1  r: 


- -fi 


+ 2K  Xn 


\(20) 


, E mat  / t \ r 

aP  -fSmain  (5^) ^ V ' V + 2K  in  ^ + 2K 

where  is  the  radios  of  the  elastic -plastic  Interface. 
SOLUTIONS  DURING  UNLOADING 


The  final  solutions  for  the  distributions  of  stresses  and 
strains  around  the  hole  are  found  by  considering  the  unloading 
phase.  During  this  phase  the  velocity  of  the  hole  surface 

t 

undergoes  a linear  deceleration  (i.e.,  u • m t)  from  u^  to  zero. 

The  derivation  of  the  equation  is  similar  to  the  loading  and 
constant  velocity  cases  discussed  earlier. 

In  the  derivation  of  the  equations  of  stress  and  strain 
distribution  for  the  constant  loading  phase  we  have  used  the  velocity 
stress,  and  strain  conditions  at  t ■■  t and  r a as  initial 
conditions;  similarly,  we  use  the  conditions  at  t ■ t2»  and  r^  • a 
as  initial  conditions  for  the  unloading  phase.  The  stress  and 
strain  distributions  at  the  end  of  the  unloading  phase,  u - 0, 
are  the  conditions  that  will  be  used  in  the  next  section  to  derive 
the  residual  stresses  and  strains. 
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For  the  elastic  region  where  r^  ^ R we  have;  (see  Appendix 


D) 


I 

- ptn'ain  (^)  - F (t) 


- pm  ain  (~)  + F (t) 


’e 


(21a 


F (t) 


- TT  F (t) 


(21b 


where 


t • 


mt 


F(t)  - mt^t  + ^ +njt^t2  ’ "”2^  ' where  t • t-t2(22) 


For  the  plastic  region  of  a ^ r^  ^ R the  equations  for  the 
stress  distribution  are 


jP  • (2K  - pm' a)  ;gn/-f|+  pm  a^n  (“^)  - ^ *^^9^ 


\ 


(R  ) 


V (23. 


■ (2K  “ P»n'a)  £n[-fj+  pm' a£n  (-^)  - £-a£i|l  + 2K 

iR  / R ) 


The  equations  of  strain  distribution  are  given  by 


3.  RESIDUAL  STRESSES:  DERIVATION  OF  EQUATIONS 


DERIVAT/.ON  OF  EQUATION  FOR  INTERFACE  PRESSURE 


Although  the  velocity  at  the  end  of  the  unloading  phase  is 
zevo,  the  final  state  of  equilibrium  between  the  rivet  and 
surrovinding  materia?  does  not  occur  until  the  materials  relax 
under  the  existing  internal  forces.  The  degree  of  relaxation 
varies  according  to  the  growth  of  the  elastic-plastic  radius 
during  loading  and  also  the  Young’s  modulus  of  each  material. 


Let  us  assume  that  in  the  final  state  of  equilibrium,  the 

interface  pressure  between  the  rivet  and  the  hole  is  P.  The 

change  of  the  radial  stress  at  r^  * a (the  rivet-hole  interface) 

is  AP  * ( cr^  -P) , where  is  given  by  Eq.  (23)  evaluated  at  r^  ■ a 

and  t ■ t^  (the  end  of  the  velocity  cycle).  The  radial 

displacement  at  r^  » a,  due  to  the  unloading  of  the  material 

surrounding  the  rivet,  assuming  incompressibility,  is  given  by 

(Ref.  8,  p.  110)  « 

(1-Hi  ) b^a 

Ad^  = . ( „ f.p)  —A (24) 


where  the  subscript  1 is  in  reference  to  the  material  surrounding 
the  rivet,  and  b>>  a. 


At  the  end  of  the  velocity  cycle,  when  the  velocity  is  zero, 
the  rivet  is  in  dynamic  equilibrium.  The  stress  deviators  that 
caused  the  plastic  deformation  in  the  rivet  during  the  velocity 
cycle  are  zero.  As  a result  of  these  conditions  the  rivet  is 
still  under  an  axial  stress  from  the  Stress  Wave  Tool  that  is 
equal  to  the  radial  stress  cr^  , the  load  from  the  surrounding 
material:  or,  in  other  words,  the  rivet  before  unloading,  is  under 
hydrostatic  pressure  a P , 
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It  can  be  shown  that  if  the  rivet  is  considered  separately 
as  a free  body  under  hydrostatic  pressure  , the  rivet  does 
not  change  in  size  when  is  removed.  Therefore,  when  the  load 
due  to  the  Stress  Wave  Riveter  is  removed,  the  radial  change  in 
the  rivet  is  due  to  the  unloading  of  the  surrounding  material, 
or  to  the  external  pressure  P.  Thus  the  radial  displacement  of 
the  rivet  during  unloading  is  given  by  (Ref.  3,  p.  322) 

‘‘^2  - - i:  (1  - (25) 

where  the  subscript  2 is  in  reference  to  the  rivet  material. 

The  condition  of  compatibility  at  the  rivet-hole  interface, 
after  complete  unloading,  is 


Ad^  - Ad2  or 


- ( 


(l+a|^)b^a 

„ .,2  27 

Ej^(b  -a  ) 


P (1  - ^i2)a 


E 


2 


(26) 


The  stresses  at  the  end  of  che  unloading  period  are  given  by 


+ Aa 


and  o-£,  * cr?  + 


A a. 


(27a) 


(27b) 


r 


— P 


r -a 
o 


k2  2 

b -a 


(27b) 


The  values  of  o ^ and  o are  given  by  Eq.  (23)  for  t * t-* 

The  circumferential  strain  at  the  end  of  phase  3 (i.e.  t • t^) 
is  found  from  Eqs.  (23b). 

The  change  in  circumferential  strain  for  the  unlcading 
phase  is  given  by 


{'■'“.V}  (^) 


and  the  change  of  radial  strain  is  given  by 


" - -a 

o . 


Thus,  the  final  residual  strain  distributions  are  given  by 


■ 'a 


(^r) 


€ + Ae 

r ^ r 


(28a) 


(28b) 


CONDITIONS  OF  UNLOADING  WITH  YIELDING 


Prager  and  Hodge  (Ref.  8,  p.  Ill)  treat  this  case,  so  we  need 
only  incorporate  their  solution.  The  criterion  for  yielding  is 
the  radial  stress  change  I (o^^)  — pl  ■ AP , where  (o^P)  is  obtained 
from  Eq.  (23a).  ^ ^ 


1') 


When  Ap  > 2K  ( 1 - — ^ ),  yielding  occurs  (Ref.  8,  p.  105). 

b 

The  value  of  ^P,  derived  from  elastic  conditions  of  unloading, 

will  provide  an  upper  boimd  value  for  the  radius  of  the  elastic- 

plastic  interface.  This  assumption  is  justified,  since  the 

value  of  Rrt  is  close  to  the  value  of  r "a. 

2 o 

The  value  of  R2  can  be  calculated  from  (Ref.  8,  p.  110) 

2 


2 in  . + -4- 
R"^  b 


2 • 


For  b>>R2  we  can  neglect  I -jg- 


Thus,  we  obtain 


R2  ' be 


-1/2  (1+2  to  I - II  ) 


(29) 


If  R2  > a then  yielding  is  occurring  during  unloading;  if  R2  < a 
then  yielding  does  not  occur  and  Eqs.  (27)  and  (28)  may  be  used. 

In  the  case  where  R2  > a we  have  the  following  equations; 


R. 


2K  (1  - *4“  ” 2 in  ^ ) a ^ r^  ^ R2 
b z 


A(J 


R.  ,2 

-2K  (1  - 


o 


R„  ^ r ^ b 
z o 


(30a) 


(30b) 


R 4. 

-21.  (1+^+2  n|^) 


A a. 


R.  .2 

•2K  -f-  (1  + -“ 


Ro  ^ r b 
Z o 


(31a) 


(31b) 


Ae, 


KR2  (1  + a) 

E r 2 
o 
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(32a) 


(32b) 


Ae  = 
r 


(1+w) 


E r 


The  final  residual  stresses  are  found  by  combining  Eq.  (23a) 
with  Eq.  (31)  and  < / 


0 * + Aa 

r r r 


‘Q 


(33  a) 


The  final  residual  strains  are  found  by  combining  Eqs.  (23b) 
and  (32) 


The  residual  circumferential  strain  is  found  by  combining 
Eqs,  (21)  and  (32).  The  following  flow  charts  show  the  sequence 
of  computations  for  calculating  the  residual  stresses  and  strains. 

FLOW  CHARTS  FOR  CALCULATIONS  OF  RESIDUAL  STRESSES  AND  STRAINS 

a)  Stresses  and  strains  at  the  end  of  the  load:Jng  phase 


Eq.(23a) 


Eq. (23b) 
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ii 


b)  Residual  stresses  and  strains  at  the  end  of  the  unloading 
phase 


P Eq.  (26) 


Eq.  (29) 


No 

R2  > a 

Yes 

r 

> R2 

0 

' 

— ] 

[IJo 

^r  * 

Q 

CD 

pEq.  (27a) 

a 

r 

> 

■ 

> _ , 

t 

SI 

[•Eq.  (28) 

m 

m 

Yes 


Eq. (30b) 
(31b) 
(33a) 


, eg  -Eq.(33b) 


(31a)“ 

(33b) 


l-Eq.  (33b) 


CONDITIONS  FOR  COMPLETE  UNLOADING  OF  THE  HOLE 

The  equations  derived  in  the  previous  sections  can  also  be 
used  for  the  case  where  the  rivet  is  removed  after  the  normal 
rivet  installation.  Such  cases  are  represented  by  selecting  the 
elastic  modulus  for  the  rivet  E2  “ 0. 

The  stress  and  strain  distributions  for  complete  unloading 
by  removal  of  the  rivet  (i.e.,  radial  stress  at  the  hole  surface 
is  zero)  are  different  from  the  distributions  with  the  rivet  in 
place.  The  differences  and  their  significance  are  discussed  in 
Section  5. 


Axial  stresses  are  calculated  from  equations  derived  in 

Appendix  E.  Axial  stresses  are  of  special  interest  in  the  cases 

where  the  rivet  is  removed  and  the  radial  stress  at  r ■ a is 

o 

zero;  then,  the  value  of  the  residual  hoop  stress  may  be  higher 

than  (jy.  This  does  not  mean,  however,  that  we  have  violated  any 
assumptions  since  we  must  also  consider  the  axial  stress  distributions. 
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COMPUTATIONS 


The  four  important  parameters  in  our  calculations  are: 

a)  Residual  radial  stress 

b)  Residual  hoop  stress 

c)  Residual  hoop  strain 

c)  Residual  radial  strain. 

The  following  cases  can  be  computed  with  data  already  in  the 
program . 


1) 


2) 


Sheet  Material 


Rivet  Material 


7075-T6 


The  programs  will  print  out  and  plot  the  radial  distance 
versus  the  residual  stresses  and  strains;  then  by  setting  E2  "■  0 
the  residual  stresses  and  strains  for  complete  unloading  can  be 
found.  The  input  necessary  for  the  running  of  the  program  is 
the  time  history  of  the  velocity  at  the  hole  surface  discussed 
in  Section  4 under  Experimental  Velocity  Measurements, 


^cedj:y(^ 

• _ .• . ‘ '»"■>“■ 

4.  EXPERIMENTAL  WORK 


INTRODUCTION 

There  are  three  separate  types  of  experimental  work  that 

have  been  carried  out  in  this  program.  The  first  consists  of 

measuring  the  force  rate  in  an  instrumented  11  foot  long 

steel  bar  (Ref.  4).  This  technique  is  equivalent  to  the 

Hopkins on ' s bar  used  in  earlier  work  with  stress  waves  (Ref.  9). 

The  second  type  of  experimental  work  consists  of  measuring  the 

' "k 

residual  strain  photoelastically  using  the  Moire  method.  The 
results  from  these  tests  are  compared  with  the  theoretical  results 
discussed  earlier.  In  the  third  type  of  experiment,  fatigue 
specimens  with  precracked  holes,  riveted  with  the  Stress  Wave 
Riveter,  were  tested  (with  precracked  unriveted  controls)  to 
show  the  effects  of  residual  stresses  in  improving  the  fatigue 
life  of  rivet  joints. 

EXPERIMENTAL  VELOCITY  MEASUREMENTS 

The  time  intervals  t^,  and  t^  were  determined  experi- 

mentally by  instrumented  bar  measurements.  A rivet  was  placed 
between  the  Stress  Wave  Riveter  and  the  11  foot  bucking  bar 
instrumented  with  strain  gauges  (see  Fig.  2).  As  the  rivet 
deforms,  the  bar  records  the  time  history  through  the  strain 
gauges  as  the  stress  wave  travels  along  the  bar.  The  following 
explanation  of  the  features  of  the  measurements  shown  in  Fig.  2 
is  based  on  our  present  understanding  of  the  forming  process. 

•k 

Grumman  Report  MMT  -379,  "Moire'  Fringe  Analysis  of  the  Strain 
FieJ.dt;  Surrounding  a Stress  Wave  Driven  Rivet,"  Engineering 
Test  Operations,  November  1977. 
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Initially,  the  radius  of  the  rivet  expands  freely  inside 
the  hole.  Since  the  cross"sectional  area  increases  during 
deformation,  the  force  transmitted  to  the  bar  increases  with 
time  because  the  plastic  modulus  is  roughly  constant.  When 
the  rivet  Impacts  the  hole  surface,  there  is  a momentary  drop  in 
force  rate  due  to  wave  reflection  in  the  rivet.  In  the  particular 
case  shown  in  Fig.  2 (A-286  rivet  in  aluminum)  the  rivet  impact 
shows  as  a drop  in  force  rate.  In  this  case,  the  acoustic 
impedance  of  the  rivet  in  the  plastic  state  is  smaller  than  the 
acoustic  impedance  of  the  aluminum.  As  soon  as  the  aluminum 
becomes  plastic,  however,  the  rate  at  which  the  force  increases 
is  restored.  (Aluminum  or  titanium  rivets  may  show  different 
behaviors . ) 

The  results  also  show  that  the  rate  of  force  increase  with 
time  is  approximately  constant,  which  indicates  a constant  rate 
of  increase  in  rivet  cross-section  if  the  plastic  modulus  is 
constant.  When  the  rivet  is  formed  there  is  no  further  increase 
in  force,  because  there  is  no  further  increase  in  area,  even 
though  the  axial  stroke  continues.  Only  when  the  force  associated 
with  forming  the  tail  of  the  rivet  is  transmitted  through  to  the 
bar  is  there  a slight  increase  in  the  indicated  force. 

The  time  measurements  from  the  typical  case  of  Fig.  2 are 
used  to  find  the  velocity  history  at  the  hole  surface  ( r » a) . 
From  previous  work  (Ref.  4)  it  is  known  that  the  acceleration 
during  the  loading  and  deceleration  during  the  unloading  phase 
can  be  approximated  by  a linear  function  of  time.  The  slope 
of  the  function  is  found  experimentally  by  measuring  the  times 
of  acceleration  (or  deceleration)  and  the  maximum  velocity.  The 
maximum  radial  velocity  is  found  from  the  following  equation: 


26 


where 


a force  rate  measured  at  the  11  foot  bar 
m 

Kj.  spring  constant  of  the  rivet  during  plastic  deformation. 

1.18  X 10^  7^ 

in. 

and 


Up  Poisson' s ratio  during  plastic  flo^l7  - 

From  Figure  2,  change  of  strain  during  t - (t.  -t  ) is 

Ae  « 100  X 10  , and  At  • 120  x 10  sec.  Thus,  a • 

sec  ’ m 


A£ 

At 


E • 7 - 19.63  x 10^  and  u - 83.2  — . 
4 r sec 


Since  during  plastic  flow  of  the  rivet,  the  velocity  u^. 
remains  essentially  constant  (Ref.  4)  an  alternate  method  to 
calculate  the  velocity  is  to  measure  the  final  expansion  of 
the  rivet  radius  so  that 


u 


r 


- ^ 

At 


For  the  case  considered  in  Fig.  2,  AR  - 0.0095  in. 
Hence, 


u 


r 


- 80 


in. 

sec 


The  two  values,  obtained  independently,  agree  well.  The 
value  of  u^  is  the  input  into  the  computer  program,  corresponding 
to  the  radial  expansion  that  determines  the  residual  stresses  and 
strains . 
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EXPERIMENTAL  MEASUREMENTS  OF  RESIDUAL  STRAINS 


Residual  strain  can  be  measured  only  through  nondestructive 
methods.  The  measurements  used  In  this  work  were  made  by  the 

t I 

Moire  technique.  To  facilitate  the  proper  use  of  the  Moire  method, 
the  riveted  plate  with  the  appropriate  grid  coating  was  protected 
during  riveting  In  either  side  by  a ^ In.  plate  that  was  removed 
after  riveting,  leaving  the  shank  of  the  rivet  Inside  the 
coated  plate  where  the  residual  strain  measurements  are  made  (see 
Fig.  3).  After  the  rivet  was  Installed  dynamically  with  the 
Stress  Wave  Riveter  head  and  tall  of  the  rivet  were  carefully 
machined  off  in  a lathe  until  the  protective  plates  dropped  off, 
leaving  the  center  plate  with  the  rivet  chank  in  place.  The 
plate  was  then  used  In  the  laboratory  to  measure  the  residual 
strains.  Having  obtained  the  residual  strains,  the  rivet  was 
carefully  drilled  out  without  touching  the  hole  surface.  The 
residual  strain  was  measured  again  (sec  footnote  pg.  25). 


PROTECTIVE  COATED  PLATE  COATED  PLATE  WITH  RIVET 

PLATES  SHANK  IN  PLACE  AFTER 

PROTECTIVE  PLATES  WITH 
HEAD  AND  TAlu  OF  RIVET 
ARE  REMOVED. 


Fig.  3 Spacially  Designed  Specimen  for  Measuring  Residual  Strain  Using  Moire'  Technique 
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We  obtained  measurements  first  with  the  rivet  In  place  and 
second  with  the  rivet  removed  which  represented  complete  un- 
loading of  the  hole.  The  difference  of  residual  strains  In 
these  two  cases  Is  substantial,  especially  near  the  hole  surface. 

The  hole  diameter,  after  the  rivet  Installation,  was  measured 
with  a magnifying  lens.  By  comparing  the  hole  size  with  the  initial 
diameter,  we  could  find  the  hole  radial  expansion,  or  amount  of 
interference.  Two  materials  were  used  for  testing;  7075-T6  and 
2024-T3.  The  specimens  were  6 x 6 x ^5  In.  and  the  hole  had  a 
nominal  diameter  of  0.230  in. 


29 


^ace<{x^i$ 

-,..-^-7-.;.-  - — . ■■ — 


5.  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

Residual  strains  were  obtained  with  7075"T6  aluminum. 

The  computed  and  experimental  value?;  are  shown  in  Figs.  4,  5, 
and  6,  tor  radial  expansions  of  0.0095,  0.004,  and  0.002  in., 
respectively.  The  plots  show  excellent  agreement  between 
theoretical  and  experimental  values  of  hoop  strain. 

The  corresponding  hoop  and  radial  stresses  from  the  theory  are 
shown  in  Figs.  7 and  8.  Both  the  hoop  and  radial  residual 
stresses  are  compressive  near  the  hole  for  the  three  cases  of 
radial  displacement. 

An  important  characteristic  of  these  results  is  that  the 
peak  positive  (tension)  hoop  stress  remains  almost  constant  at 

— ■ 0.4.  An  excessive  level  of  would  lead  to  accelerated 

O'  0 

y 

failure  when  tension  loads  are  applied.  The  position  of  the 
peak,  however,  increases  radially,  away  from  the  hole,  as  the 
same  specimens  with  the  rivet  removed  are  shown  in  Figs.  9,  10, 
and  11.  Comparison  of  experimental  and  theoretical  results 
for  one  case  of  an  open  hole  (i.e.,  AR  - 0.0045  in.)  shows  good 
agreement  (see  Fig.  12). 

In  another  set  of  calculations  and  experimental  measurements, 
2024-T3  material  was  used  with  A-286  rivets.  The  results  are 
shown  in  Figs.  13,  14,  and  15,  (i.e.,  AR  « 0.011  in.).  The  case 
of  2024-T3  specimens  with  the  rivet  removed  is  shown  in  Figs. 

16,  17,  and  18. 


HOOP  STRAIN 


HOOP 


MATERIAL  7076.T6 
RIVET  1/4  IN.  A-286 
AR  - 0,002  IN. 
4^EXPERIMCNT 
— THEORY 


(!!■ 


r 


DIMENSIONLESS  RADIUS  Rn/A 
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Fig.  10  Residual  Radial  Stress  vs  Radial  Distance 


Fig.  15  Residual 


/4  IN.  A286 


6.  CONSTANT  AMPLITUDE  FATIGUE  TESTS  OF  SPECIMENS  WITH  RESIDUAL 


STRESSES 

When  residual  stresses  are  combined  with  applied  stresses 

due  to  a uniaxially  applied  load  in  a specimen,  they  directly 

affect  the  fatigue  life  of  the  specimen.  To  illustrate  this 

effect  specimens  made  of  2024~T87  material  were  used.  In  the 

first  series  of  tests  specimens  whose  geometry  is  shown  in  Fig.  19 

were  made  each  with  three  open  holes  of  0.188  in.  diameter. 

They  were  tested  at  constant  amplitude  and  the  results,  along 

*ic 

with  previous  results,  are  shown  in  Fig.  20. 


HOLE  DIAM. -3/16IN. 


Fig.  20  Fatigui  Lift  of  202<i  T81  Aluminum  Alloy  Optn  Holt  tnd  Counttnunk  Optn  Holt  ConUol  Sptcimeni 


Note:  Grumman  Report  MPLR  45-76-23,  "Constant  Amplitude  Fatigue 
Test  Results  of  Cracked  and  Uncracked  2024-T81  Aluminum 
Alloy  Containing  "Unloaded"  Stress  Wave  Driven  Rivets," 
Materials  and  Process  Department,  19  May  1977. 
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In  the  second  series  of  tests  the  specimens  were  drilled 
with  the  same  holes  and  riveted  with  the  Stress  Wave  Riveter 
using  three  different  types  of  rivets;  2117-T4,  A-286  (GR-501W-6), 
and  Beta  C titanium.  The  specimens  were  again  tested  at  constant 
amplitude  at  two  stress  levels.  30  ksi  and  20  ksi.  The  results 
are  shown  in  Figs.  21,  22,  and  23.  The  radial  displacement  was 
0.005  and  0.006  in.  in  each  case. 

The  important  conclusion  from  this  limited  number  of  speci- 
mens is  that  regardless  of  rivet  material,  the  improvement  in 
fatigue  life,  when  compared  to  open  hole  results,  is  about  ten 
times  that  of  the  open  hole  specimens.  More  extensive  testing 
is  required  to  arrive  at  more  definitive  conclusions.  It  is 
also  important  to  find  some  theoretical  mechanism  to  predict 
life  once  the  initial  residual  stresses  and  the  load  are  known. 

A modified  Goodman' s diagram  may  be  suitable. 
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STRESS  (S),  KSI 


Fig.  22  Fatigue  Life  of  2024-T81  Aluminum  Alloy  Speoiment  Containing  Unloaded  3/16  In.  Dia  A-286  Steel  Rivets 
Installed  with  Stress  Wave  Riveter. 


O UNCRACKED  SPECIMENS 


OO  00 


(2)  NFO 


FATIGUE  LIFt-'  (N),  CYCLES 

Fig.  23  Fatigue  Life  of  2024-T81  Aluminum  Alloy  Specimens  Containing  Unioaded  3/16  In.  Dia  Beta  C Titanium  Rivets 
Installed  with  Stress  Wave  Riveter 


7.  EFFECT  OF  PRECRACKED  HOLES  ON  FATIGUE  LIFE  OF  DYNAMICALLY 
RIVETED  SPECIMENS 

A radial  crack  in  a hole  surface,  where  stress  concentration 
is  high,  provides  the  most  favorable  conditions  for  failure.  To 
show  the  beneficial  effect  of  dynamic  riveting  using  the  Stress 
Wave  Riveter,  we  tested  specimens  with  precracked  holes. 

The  specimens  were  of  the  same  design  as  described  in 
Section  6.  Pilot  holes  of  ^ in.  were  first  drilled  and  the 
specimens  were  subjected  to  constant  amplitude  fatigue  until  a 
radial  crack  was  developed  in  a preassigned  notched  area.  The 
holes  were  then  drilled  and  countersunk  to  il  in.  diameter, 
with  the  crack  in  each  hole  of  0.05  in.  measured  to  within 
0.005  in.  and  extended  through  the  thickness  of  the  countersunk 
part.  (Each  specimen  consisted  of  two  parts,  i in,  thick  each 
made  of  2024-T81  material.) 

The  Stress  Wave  Riveter  was  used  to  install  rivets  made  of 
(a)  modified  Beta  C Titanitom  rivets  (b)  2117-T4  aluminum  rivets, 
and  (c)  A-286  rivets.  The  nominal  radial  expansion  was  0.006 
in.  in  each  case. 

The  specimens  were  tested  at  constant  amplitude  load 
cycling.  The  results  are  shown  in  Fig,  24.  The  specimens 
riveted  with  A"286  show  the  best  improt'ement  when  compared  to  the 
fatig;ue  life  of  open  hole  precracked  specimens.  It  is  apparent, 
however,  that  all  the  riveted  specimens,  with  compressive  residual 


Note:  Grumman  Report  MPLR-45- 77*135 , "Open  Hole  Precracked 
Specimens  Tested  at  Constant  Amplitude  Load  Cycling," 
Materials  and  Process  Department,  October  1977. 
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stresses,  improved  fatigue  life  and  otherwise  slowed  the  crack 
growth.  The  fatigue  life  of  each  riveted  specimen  was  approx- 
imately ten  times  greater  than  the  open  hole  precracked  specimens. 


OhK  o 


MATERIAL  2024-T8i  SPECIMENS 
RIVETED  WITH  STRESS  WAVE 
RIVETER  AT  0.006  IN.  INTERFERENCE 
EACH  HOLE  PRECRACKED  PRIOR 
TO  RIVETING  WITH  0.06  IN.  CRACK 


! 10 


X OPEN  HOLE  PRECRACKED 


. PRECRACKED  & RIVETED  WITH 
BETA  C TITANIUM  RIVETS 


I PRECRACKED  & RIVETED  WITH 
'2177-T4AL  RIVETS 


Q PRECRACKED  & RIVETED  WITH 
A-286  RIVETS 


FATIGUE  LIFE  (N),  CYCLES 

Fig.  24  Constant  Amplitude  Load  Cycling  Fatigue  Tests  cf  Specimens  Riveted  with  Stress  Wave  Riveter 


8.  DISCUSSION  OF  RESULTS 


The  plane  strain  assumption  made  in  the  analysis  has  given 
results  that  agree  well  with  tests.  It  must  be  remembered, 

I 

however,  that  the  specimens  used  with  the  Moire  technique 
were  ^ in.  thick  and  were  protected  during  riveting  on  either 
side  by  in.  plates.  Metallographic  analysis  shows  that  near 
the  surface  plane  stress  conditions  prevail.  These  conditions 
were  removed  with  the  protective  plates.  The  plane  strain 
conditions  therefore  are  not  valid  in  thin  sheets,  less  than 
0.06  in.  thickness,  where  the  present  analysis  may  not  be 
applicable. 

The  plane  strain  condition,  of  course,  is  accompanied  by 
axial  stresses  as  shown  in  Appendix  E.  With  the  rivet  in  place, 
triaxial  stresses  exist  near  the  hole  surface.  When  the  rivet 
is  removed,  however,  we  have  only  hoop  and  axial  residual  stresses. 
These  conditions  can  be  used  with  a modified  Goodman' s diagram 
to  predict  fatigue  life. 

The  incompressibility  assumption  was  essential  in  the 
development  of  the  theoretical  analysis.  Prager  and  Hodge 
(Ref.  8)  argue  that,  "the  assumption,  which  greatly  facilitiates 
the  mathematical  work,  furnishes  good  results  if  the  corrections 
2u  and  2(l-u)  ara  applied  to  the  axial  stress  and  the  radial 
displacement,  respectively".  They  arrived  at  these  factors  by 
comparing  theoretical  solutions  of  compressible  and  incompressible 
materials  in  static  load?  ig.  Since  we  are  comparing  the  in- 
compressible solution  with  experimental  results  in  dynamic  loading 
we  believe  that  any  discrepancies  are  within  the  errors  of 
computation  and  therefore  no  correction  factors  are  needed. 
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The  radial  and  hoop  stress  solutions  between  compressible  and 
incompressible  materials  should  agree  well  according  to  Prager 
and  Hodge. 

I 

The  Moire  technique  used  to  measure  residual  strains  proved 
successful  enough  to  recommend  it  for  practical  applications. 
Despite  early  difficulties  with  fringes  disappearing  near  the 
hole,  the  final  quality  of  the  fringe  reproduction  is  excellent. 
Our  Innovative  technique  of  protective  plates  during  riveting 
was  Imperative  for  the  protection  of  thin  film  near  the  hole. 

In  conclusion,  the  results  of  this  work  show  that  there 
is  a triaxial  residual  stress  field  around  the  hole  that  is 
induced  by  dynamically  expanding  the  rivet  in  the  hole.  Further- 
more, the  radial  extent  of  cold-work  is  at  least  as  much  as  that 
provided  by  other  cold-work  methods  such  as  the  sleeve  and  the 
split  sleeve  techniques.  In  addition,  it  is  shown  that  the 
residual  stresses  do  improve  fatigue  life.  Comparison  of 
fatigue  data  shows  that  for  constant  amplitude  cyclic  loading 
and  0.006  in.  radial  expansion  of  cold-work  we  obtain  an 
increase  of  10  times  the  fatigue  life  of  specimens  without 
the  cold-work.  What  is  more  important,  however,  is  that  residual 
stresses  induced  around  a precracked  hole  by  the  dynamically 
expanding  rivet  slow  down  the  crack  propagation  and  improve 
fatigue  life. 
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APPENDIX  A 

DERIVATION  OF  STRAIN  RELATIONS 
DEFINITION  OF  STRAINS 


Radial  Strain: 

Hoop  Strain: 
Volumetric  change: 
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RELATIONS  WITH  MEAN  NORMAL  STRESS 


From  Hook's  law  we  have: 
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where 
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Adding  the  three  stress  equations  we  obtain 


(CTr+a^+Og)  - 3 (l+u)  ^ • 


Since 
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we  have 
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From  the  definition  of  the  stress  deviator  we  have 

. "r+V“e  ., 

S « 2 • 3ke 

From  the  definition  of  the  principal  stresses,  we  have 


(A-4) 


• S + , cTg  - S + Sg  and  - S + 


A-2 


Adding  the  three  equations  we  obtain 


(cJj.  + cTg  + a^)  - 3S  + (S^  + Sq  + ). 

From  Eq.  (A-3),  however,  we  have  a + + o » 3S. 


Thus,  we  conclude  with 
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RELATIONS  BETWEEN  STRESS  AND  STRAIN  DEVIATORS 
From  Hook’ s law  we  have 


(1+n) 
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or 
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Substituting  for  S from  Eq.  (A“4)  in  the  right  hand  side 
and  k ■ obtain 

o ” S » (€  - e) 

r (1-Kl)  ' r ^ 

Introducing  incompressibility,  where  e • 0,  we  obtain 
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APPENDIX  B 


DERIVATION  OF  EQUATIONS  OF  STRESS  AND  STRAIN  DISTRIBUTION  DURING 
ELASTIC  LOADING 


It  is  possible  to  Integrate  each  of  Eqs.  (13)  separately 
and  evaluate  the  constants  from  the  boundary  conditions. 

c^(t) 

Integrating  the  second  equation  in  (13)  we  obtain  u ■ , 

and  from  the  boundary  conditions  of  the  hole  surface  u ■ u^  - mt 
at  r^  “ a,  we  obtain  a general  solution  for  the  velocity 
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mt»  a 
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Integrating  the  compatibility  equation  (Eq.  (2))  with  respect 

to  r using  Eq.  (3)  with  e - 0,  we  also  obtain 
o z 

C2(t) 

* 2 

^o 

Since  at  r^  » a,  u - mt  we  obtain  ^^(a)  » “ from  Eq.  (4). 

2 mt^ 

Hence  C2(t)  • ® C2(t)  ^ 2~^>  since  at  t ■ 0, 

« 0 and  €2(0)  - 0. 


The  equation  for  is,  therefore,  given  by 


mt^a 
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and  employing  the  plane  strain  condition  where  e 


■-€g,  we 


obtain 
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B-1 


The  first  equation  of  (13)  can  now  be  reduced  to  a simpler 
form.  Using  Eq.  (lOa),  we  can  write 

■^o  r 2 

o 

Substituting  in  Eq.  (10)  for  and  we  obtain 


r.1  ^2 

e . . E mt  a 

^r  2 

^ 2v  ^ 


Multiplying  by  r^  and  differentiating  we  obtain 


Thus, 


?r  <Sr'o  > - 0 


or 
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or  r 

o o 

Integrating  Eq.  (B-3)  we  have 


(B-3) 


S - - fma  ^n  r^  + C2(t) 

Where  C2(t)  is  the  constant  of  integT:atlon. 

We  finally  find  an  expression  for  S,  using  the  boundary  condition 
of  S — 0 at  all  times  when  r^  ■ b where  b » a. 

S - pma  in  ~ (B-4) 

^o 
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APPENDIX  C 

SOLUTION  OF  CONSTANT  VELOCITY  REGION 

Tests  show  (see  Fig.  2)  that  when  the  rivet  impacts  the 

hole  surface,  there  is  a momentary  reversal  in  velocity  change, 

from  acceleration  to  deceleration.  When  acceleration  begins 

again,  it  takes  time  t reach  a constant  velocity.  Then  the 

velocity  remains  constant  for  the  time  interval  (t2  " t^). 

During  this  interval  the  acceleration  is  zero  and  the  velocity 

at  r * a is  u » u . 
o o 


^ IT 

Cj.  + ~ " 0 (incompressibility) 

o ^o 

differentiating  w.r.t.  time  we  obtain 
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From  Eq.  (B"l)  we  know  that, 
mat 
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differentiating  w.r.t.  r^  we  obtain 
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and  dividing  by  r^  we  have 
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Substituting  in  Eq.  (10)  we  obtain 


and 
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The  equilibrium  equation  is  the  same  as  Eq.  (B“3)  with 
m ■>  0,  i.e.,  (constant  velocity)  that  reduces  to 


0 : or,  S ■ C^(t) 


where  is  the  constant  of  integration.  We  evaluate  from 
Eq.  (B-4)  which  applies  at  t « t^.  Thus,  S « fma  ^n  (~)  which 
indicates  that  it  is  constant  with  respect  to  time.  ° 


The  elastic  solution  during  the  constant  velocity  interval 
is  found  by  substituting  for  S^,  and  S in  the  equations  « S + 
and  CTg  - S + Sgi 
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When  the  material  becones  plastic  during  the  loading  phase, 
the  plastic  flow  will  continue  into  the  constant  velocity  phase. 

At  t * t^  the  radius  of  the  elastic -plastic  interface  R is  given 
by  combining  Eqs.  (10),  and  (B-2),  substituting  r^  ■ R, 
and  then  solving  for  R,  when  - K.  The  radius  of  the  elastic- 
plastic  interface  during  the  constant  velocity  period  E|,  is  found 
in  similar  fashion  from  Eqs.  (10),  and  (C-2),  for  « -K. 

For  r < R,  the  material  is  plastic,  for  r > R.  the  material 
o o i 

is  elastic.  The  stresses  and  strains  in  the  plastic  region 
are  found  as  in  the  case  of  increasing  velocity.  Since  we 
assume  incompressibility  during  elastic  and  plastic  cases  the 
equation  for  velocity  remains  unaffected. 

Hence , 


The  stresses  are  given  by  a » S - K and  On  “ S + K with 


+ Sj,  Sq  + Sg""  ■ K*"  and  cr^  “ » - 2K.  Substituting  into 

equilibrium  equation  (for  m «•  0)  we  obtain 
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Integrating  we  obtain 


S * 2Ki!nr  + f(t) 
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At  r^  ■ Rj^  both  the  elastic  and  plastic  solutions  hold. 


Thus, 
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Thus,  the  stresses  for  r^  < Rj^  are  given  by 
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The  stresses  for  > Rj^  are  given  by  Eq.  (C-4). 

where  R^^  « elastic-plastic  boundary  during  constant  velocity  phase, 
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where 
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Substituting  in  Eq.  (10),  we  obtain 
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Substituting  in  the  equilibrium  equations  we  obtain  from  Eq. 
(B-3) 
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^r  " r 
o o 


where  m is  the  deceleration  rate. 
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Thus,  the  elastic  solution  is  given  by 
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for  t2  ^ t ^ t^. 
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The  plasticity  condition,  S^.  “-K  and  Sg  • K will  continue 
into  the  unloading  phase.  Substituting  in  the  equilibrium 
equation  we  obtain 

^ 2j^  _ pm  a 

Sr  " r r 
o o o 

1 

After  integration  we  have  S » (2k  - pm  a)  + C^*  Since  at 

the  radius  of  the  elastic-plastic  interface  r'  , both  the  elastic 
and  plastic  solutions  hold,  we  obtain 


(2K  - pm  a)  inR  + C,.  - K « pm  a in  (“) 
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Hence, 
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Thus,  for  r^  < R , the  solution  is  given  by 
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The  value  of  R is  found  from  S - -K  - - 
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The  stresses  for  r > R are  given  by  Eq.  (D“5).  At  the  | 

o II  I 

f the  unloading  phase,  t to.  int_  •-m  t , and  u ■»  0.  ' 1 


APPENDIX  E 


AXIAL  STRESSES 


With  plane  strain  conditions,  we  have  S • 0.  Using  Eq.  (7) 

z 

we  obtain  » S.  The  final  conditions  given  in  Appendix  D 
provide  the  following  equation  for  S. 


Elastic  case 
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Plastic  condition  a^«S»a^+K  (E-2) 

z r 

The  axial  stress  for  unloading  with  yielding  is  taken  from  Ref.  8 


Ac  « Aa  ••2k 
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Thus,  the  final  axial  stress  (o)  • aj?  -f  Aa^:  or,  combining 

z z z 

Eqs.  (E-2)  and  (E“3),  we  obtain 
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When  the  rivet  is  removed,  (a)  • 0 at  r • a. 
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The  axial  stress  for  unloading  without  yielding  is  given  by 

A a -(P  - aP) 

b -a 
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or,  using  Eq.  (27),  we  have 
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Separate  calculations  are  necessary  to  obtain  axial  stress 
distributions,  since  they  are  not  included  in  the  computer  program 
discussed  earlier. 
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times  higher  than  the  fatigue  life  of  specimens  with  open  holes.  In  the 
second  series,  precracked  specimens  with  a 0.05  in.  crack  in  each  hole  were 
riveted  with  the  StressWave  Riveter.  The  cracks  were  oriented  so  that  the 
fatigue  loading  was  normal  to  the  crack  path.  The  fatigue  life  of  these 
specimens  was  also  approximately  ten  times  higher  than  the  fatigue  life  of 
those  specimens  with  precracked  holes  with  rivets. 


